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We derive Ginzburg-Landau action by systematically integrating out electronic degrees of freedom 
in the framework of the Keldysh nonlinear cr-model of disordered superconductors. The resulting 
Ginzburg-Landau functional contains a nonlocal A-dependent contribution to the diffusion constant, 
which leads, for example, to Maki- Thompson corrections. It also exhibits an anomalous Gor'kov- 
Eliashberg coupling between A and the scalar potential, as well as a peculiar nonlocal nonlinear 
term. The action is gauge invariant and satisfies the fluctuation dissipation theorem. It may be 
employed, e.g., for calculation of higher moments of the current fluctuations. 
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I. INTRODUCTION 



Time dependent Ginzburg-Landau (TDGL) theory 
has received a lot of attention and was a subject of 
controversy over many years 

Gor'kov and Eliashberg^ (GE) were probably among the 
first who realized that the thermodynamic Ginzburg- 
Landau equation may be generalized for the time- 
dependent phenomena in the case of gapless supercon- 
ductivity (see also earlier publicationsii^) . The latter oc- 
curs either in the presence of magnetic impurities, or in 
the fluctuating regime at T > T^. Notably GE equa- 
tion contained an anomalous nonlocal coupling between 
the order parameter A and the scalar potential: the fact 
that was frequently overlooked in many subsequent treat- 
ments. 

Extension of the TDGL theory to a gapped phase 
turns out to be a very demanding problem. As noted 
by Gor'kov and Eliashberg, the difficulty stems from the 
singularity of the BCS density of states at the gap edge. 
The latter leads to a slowly decaying oscillatory response 
at frequency 2A/h in the time domain. As a result, the 
expansion in powers of the small parameter A/Tc fails. 
In principle, it may be augmented by an expansion in 
A/{huj), in case the external fields are high-frequency 
ones. To describe low-frequency responses in the gapped 
phase, one needs a time nonlocal version of the TDGL 
theory. The analysis is greatly simplified in the presence 
of a pair-breaking mechanism, such as magnetic impuri- 
ties or energy relaxation. Such a mechanism may elimi- 
nate singularity in the density of states, leading to gapless 
phase in the presence of finite A. Under these conditions, 
an expansion in powers of At^/H and ujt^ is justified and 
thus a time-local TDGL equation may be derived (here 

is the pair-breaking time). In the present work, we re- 
strict ourselves to the fluctuating regime T > Tc, where 
the spectrum is gapless automatically and there is no 
need in an explicit pair-breaking mechanism. 

Soon after the GE work, Aslamazov-Larkin (AL)i^ 
and Maki-Thompson (MT)iiii^ corrections to conduc- 
tivity of fluctuating superconductors were discovered in 
the diagrammatic linear response framework. While AL 
term had naturally followed from TDGL theory (see, e.g.. 



books22i2ii^), MT phenomena were seemingly absent in 
TDGL formalism. Based on the work^^ it was proposed^ 
that in order to include MT term into the set of TDGL 
equations, one has to substitute the renormalized con- 
i.2.3.4.5.6.7 .8.9.io.ii.i2."i3.i4.i5 ductivity (T ^ (T + cr^^"^ in the expression for the current, 

supplementing TDGL equation. While leading to the 
correct static average current (by construction) , this way 
of handling the problem fails to satisfy the fluctuation- 
dissipation theorem (FDT). Indeed, it does not provide 
any prescription for calculating higher moments of the 
current (even in equilibrium). Another drawback of the 
approach of Ref. Q is that it fails to incorporate a pecu- 
liar frequency dependence of MT phenomena, stemming 
from the time nonlocality of MT terms. The procedure 
introduced phenomenologically in Ref. [§] was latter el- 
egantly derived in Ref. [l9| using nonequilibrium Green 
functions technique. Let us also mention few other works 
where a combined set of TDGL and kinetic equations was 
suggested i^^d'^ An imaginary-time action of fluctuating 
superconductors was discussed in Ref. 



15| 



In the present publication, we derive a set of coupled 
stochastic TDGL and Maxwell equations, which are suit- 
able for calculation of both average current and its higher 
moments. This set of equations is an immediate conse- 
quence of the effective Keldysh action written in terms of 
the fluctuating order parameter and electromagnetic po- 
tentials. Technically, we employ the nonlinear cr-model in 
the Keldysh representatio n^^i^^ to perform disorder av- 
eraging. We then systematically integrate out the elec- 
tronic degrees of freedom, neglecting Anderson localiza- 
tion effects. The resulting effective action, written in 
terms of the order parameter and electromagnetic po- 
tentials, naturally and unmistakably contains both MT 
terms and anomalous GE coupling between the order pa- 
rameter and electric field. 

We restrict ourselves with the fluctuating regime T > 
Tc only, leaving the case T < Tc (and magnetic impuri- 
ties) for future studies. As always, the Ginzburg-Landau 
treatment requires the condition 



T-Tc^Tc, 



(1) 



which is central to our consideration. We also assume 
that both the order parameter and the electromagnetic 
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fields v ary on the spatial scale which is much larger than 
^0 = yD/Tc (here D is the diffusion constant) and the 
time scale which is much slower than 1 /T^ (hereafter we 
adopt units, where h = c = 1). Moreover, we shall rely 
on the fact that the electronic system is always in a lo- 
cal thermal equilibrium. This in turn implies that the 
external fields are not too large. More precisely, the elec- 
tric field E is such that ei?^o ^ 7c, while the magnetic 
field H is restricted by the condition eiJ^o l/^o- No- 
tice that these conditions do not restrict our treatment 
to the linear response regime. Nonlinear phenomena may 
be included, as long as a characteristic scale of nonlinear 
effects satisfies the inequalities given above. 

The restrictions on spatial and temporal scales of the 
external fields along with the fact that electrons are in 
local equilibrium considerably simplify the theory. In 
particular, most of the terms in the effective action ac- 
quire a local form in space and time. Nevertheless, the 
effective theory does not take a completely local form. 
The diffusion constant obtains a A-dependent contribu- 
tion, with essentially nonlocal coupling to the order pa- 
rameter. If averaged over the fiuctuations of the order 
parameter, this nonlocal term yields MT correction to 
conductivity. We note, however, that an average current 
is not the only manifestation of the nonlocal term. The 
latter also contributes to the current noise as well as to 
its higher moments. Another nonlocal effect in the effec- 
tive action is the way the order parameter interacts with 
the time-dependent electric field. This is the anomalous 
GE term. There is one special gauge (A^-gauge), where 
an anomalous term takes an especially simple form. In 
what follows, we shall explain the M^-gauge and perform 
all the calculations in it. The resulting action may be 
then transformed back into an arbitrary gauge. 

The use of the Keldysh formalism is important in sev- 
eral respects. First, it allows to augment the replica 
trick to perform the quenched disorder averaging proce- 
dure. Second and more important, it is the only consis- 
tent way to derive real-time dynamics. The use of the 
imaginary-time formalism, although possible, requires 
performing the analytical continuation procedure. The 
latter is known to be exceedingly demanding for MT as 
well as time nonlocal nonlinear terms. Working directly 
in real time allows to make all the expressions physically 
transparent, unobscured by the peculiarities of the ana- 
lytical continuation. Finally, the Keldysh formalism nat- 
urally allows to extend the treatment to the situations, 
where the assumption of local equilibrium is not appli- 
cable. Although not considered in the present work, a 
treatment of a nonequilibrium fluctuating superconduc- 
tivity is a subject of great interest. 

The rest of the paper is organized as follows. In the 
next section, we present our main results in the form 
of the set of coupled stochastic equations for the order 
parameter and electromagnetic potentials. In Sec. IIIII 
we introduce the basic elements of the Keldysh nonlinear 
CT-modcl and explain the way the effective action is de- 
rived by integrating out diffuson and Cooperon degrees 



of freedom. Technical details of this procedure are dele- 
gated to a number of appendixes. Finally, in Sec. llVi we 
summarize our findings and briefly discuss their possible 
applications. 



II. SET OF STOCHASTIC EQUATIONS 

The most compact way to present our results is in the 
form of the effective Keldysh action which is a functional 
of the fluctuating order parameter and electromagnetic 
potentials. Since it requires introducing some notations, 
we postpone discussion of the action until Sec. IIIII Here, 
we present an equivalent way to display the same infor- 
mation using the set of stochastic TDGL and Maxwell 
equations. 

In presence of the scalar ^(r, t) and vector A(r, t) po- 
tentials the complex order parameter A(r,i:) obeys the 
following TDGL equation: 

{dt - 2ie dtK) A = [d ( V - 2ieAf - r'^ 1 A + Ca , (2) 

here D is the diffusion constant, e is the electron charge, 
and 



tgl = 



(3) 



is the Ginzburg-Landau relaxation time. The field 
/C(r,t) satisfy the following equation: 

{dt ~ Oy^) /C(r, t) = $(r, i) - D divA(r, t) . (4) 

The complex Gaussian noise ^A(r, ^) has the correlator 

(eA(r, t) a(r', t')) = <5r_r' St-t' , (5) 

where v is the density of states. Unlike TDGL equation 
frequently found in the literature ji'^ii^'^ii^ the Ihs. of 
Eq. ([2]) contains GE anomalous term^ dtlC{r,t) instead 
of the scalar potential $(r, i). The two coincide in the 
limit of spatially uniform potentials, cf. Eq. In a 
generic case, they are rather distinct and /C(r, t) is a non- 
local functional of the scalar and the longitudinal vector 
potentials. The standard motivation behind writing the 
scalar potential ^{r,t) in the Ihs of TDGL is the gauge 
invariance. Notice, however, that a local gauge transfor- 
mation 



A ^ Ae-2«x. ^^<i>-dtx; 
A ^ A-Vx; /C^/C-x 



(6) 



leaves Eq. unchanged and therefore this form of 
TDGL equation is perfectly gauge invariant. The last 
expression in Eq. © is an immediate consequence of 
Eq. ([4]) and the rules of transformation for <I'(r,t) and 
A(r,t). 

We have suppressed the nonlinear terms in Eq. 
since they are of lesser importance for T > Tc. A detailed 
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discussion of the nonlinear terms is presented in section 
Sec. IIIIDI We note, however, that in addition to the 
conventional |ApA local term, there is other essentially 
nonlocal and time-dependent nonlinear term in TDGL 
equation. 

TDGL equation ([2]) takes an especially simple form 
in the ]C-gauge, which is obtained by choosing x(r, t) = 
/C(r,t) in Eq. In other words, the gauge is specified 
by the relation 



<i>r - Z) divAr = . 



(7) 



where — ^ - dtX and = A - VIC. In such a 
gauge, the anomalous term in the Ihs of TDGL is absent 
and the latter obtains the form 



dtA 



D{V- 2ieA^f - T^l 



+ a , (8) 



where A^; = Ae^^*'^'^. Employing Eq. ([7]) along with 
the expression for the electric field E = dtAjc ~ V^ic, 
one finds for the vector potential A^(r,t) in the rhs of 
TDGL equation (0, 



AK(r,i) = Ai(r,i) + A||K;(r,t). 



(9) 



Here, A^ is the gauge invariant transverse part of the 
vector potential and the longitudinal part in the /C-gauge 
is given by 



A||K(r,t) = Jdr'dt'Vl-^'-E\\{r',t')., 



(10) 



where E|| is the longitudinal part of the electric field. The 

kernel V^l^:' - e{t - t') 
the diffusion operator 



kernel 'D^^, ^ 6{t — t') is the retarded Green function of 



{dt-DV^)Vll, ^St-t'Sr-r'. 



(11) 



In addition to the equation for the order parameter, the 
complete theory must provide two material equations for 
the current j(r,t) and charge p(r,t) densities. The first 
of these equations is the continuity relation: 



div j + dtp^O. 



(12) 



As for the second one, we found the following expression 
for the current density: 



/IT P 1J I ) 
dt' [DSt-r + SD^^Z'] P^^E(r,t')-Vp(r,i')] +^^Im[A;^(r,t)(V-2*eA^)AK(r,t)]+^j(r,t). (13) 

The nonlocal part of the diffusion coefficient is the functional of the order parameter (as well as the electromagnetic 
potentials) and is given by the expression 



5D^Z,[A^] = ^ / dr'dr" €^,'^1 (r',r) A^(r",r) C:^/, , 



(14) 



with T = (< + 1')/2 . The retarded C^'/"t/ ^ 9(t — t') and advanced C^fj/ ^ 6{t' — t) Cooperon propagators are Green 
functions of the following equations: 



5, 



<^K.{v,T+) + ie$/c(r,T-) - [V - zeAK;(r,T+) - zeAK:(r,T_ 



(15a) 



-dt + ie^K.(j, r+) - ie^tcij, r_) - D [V - ieAfciv, r+) - ieA/c(r, T-)f C^'j = Sr-v'^t-t' , 



(15b) 



where r± = t ± t/2. Note that the MT term obeys causality, since 5D 



MT 



9{t — t'), and gauge invariant in view of 



Eq. Being averaged over the fluctuations of the order parameter {5D'^^'^)a, it leads to the (frequency-dependent) 
Maki-Thompson correction to the conductivity. Equation (fT3)) is more general, however, as it allows to calculate the 
higher moments of the current as well. The current fluctuations are induced by the stochastic term in the TDGL 
equation as well as by the current noise ^j(r, t) given by the Gaussian vector process with the correlator: 



(Cf (r, t) (r', t')) = Te^i. {2D5,^, + SD^^Z' + ^^'A-,) 



(16) 



guaranteeing validity of FDT, here a,(3 — x,y,z. Equations (|12p . and must be supplemented by Maxwell 
equations for the electromagnetic fields. In the next section, we show how these results may be derived from the 
microscopic model. 



III. KELDYSH SIGMA MODEL FORMALISM 



real time. The formalism considers the evolution along 



A. Notations and the cr-model action 



We employ Keldysh technique^ which allows to go 
beyond the linear response and is formulated directly in 
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the closed contour in the time direction. It thus deals 
with the two "replica" of each field, one encoding the 
evolution in the forward and another in the backward 
time direction. It is convenient to introduce half-sum 
and half-difference of these fields to which we shall refer 
as classical and quantum components correspondingly.— 
As a result, all the fields acquire the vector structure, e.g., 
the scalar potential $ = (<I>^', the vector potential 
A = (A'^',A'') and the complex order parameter A ~ 
(A"^', A*). It is also convenient to introduce 4x4 matrix 
notations for these fields in the space which is a direct 
product of Keldysh and Nambu spaces 

$ = [$'='CTo + «>'?crj0To, A = [A'^'ao + AV^l^ro, (17a) 



A = [A^'ao + A V,] ® r+ - [A*='ao + A* V,] ® t_ , (17b) 

here <7a and Tq arc sets of Pauli matrices in Keldysh 
and Nambu spaces correspondingly (a = 0, x, y, z) and 

T± = {Tx ± iTy)/2. 

Our starting point is the nonlinear cr-model )^'^'^^ which 
systematically takes care of the elastic disorder averag- 
ing. In the framework of this formalism, the electron 
dynamics is described by the field Q{r;t,t') which is a 
matrix in the 4x4 Keldysh-Nambu space as well as an 
infinite matrix (integral kernel) with respect to its two 
time indices. For a short-range correlated disorder (the 
only case considered here) the Q matrix is a local func- 
tion of the spacial variable r. The Q matrix obeys the 
local nonlinear constraint 



Q'(r)-1: 



(18) 



where Q ^ is understood as the matrix multiplication in 
4 X 4 as well as in the time space and the rhs is the unit 
operator in this space. 

It is very convenient^^i^ii^ to single out the gauge de- 
gree of freedom of the Q-matrix field by expressing it in 
the following form: 

Q(r; t, t') = e-2K(r,t) Q^^^. ^/^ ^-^eK(,r,t')S (^g) 

Here A^(r, t) is a scalar gauge field in the matrix represen- 
tation analogous to Eq. p7ap . S = ctq (8)Tz, and Q/c is the 
new Q-matrix field free from the gauge ambiguity. Ob- 
viously, the field Q/c a-lso satisfy the nonlinear constraint 
= 1. In what follows, we shall use the freedom of 
choosing the gauge field AC to adjust a saddle point on 
the Q-manifold, Eq. (fT5|) . according to local scalar and 
vector potentials. Therefore, the IC field should be un- 
derstood as a certain functional of the electromagnetic 
potentials $ and A which fixes a special gauge. 

The Keldysh nonlinear a-model, we employ here, was 
formulated for normal metals by Kamenev and An- 
drce\i^ and extended for superconductors by Feigelman 
et al^ Its action takes the following form: 



2 

= Tr [AkTAk] , 5$ = ^ Tr [$kT$k] , 

(20b) 

5,= ^Tr [D{drQKf~mQK + '^^e^KQK + '^^^icQK] , 

(20c) 

here A is superconductive coupling constant, T — CTa, (Stq, 
and the covariant spatial derivative is defined according 
to 

drQic - VrQK - ie[SAK, Ok] • (21) 
The subscript K, denotes gauge transformed fields 

= $ - dtIC, Aic = A- VIC, (22a) 



AK:{r,t) =e-*'=^^('''*)A(r,t) e 



ielC{r,t)'E 



(22b) 



The trace operation in Eq. ([^D|) Tr[. . .] implies integra- 
tion over the space and time indices as well as matrix 
trace in the 4x4 Keldysh-Nambu space. The action 
written above should be supplemented by the standard 
Maxwell term Sm = Tr [ETE + HTH] /IGn. 

Our eventual goal is to integrate out fluctuations of 
the electronic degrees of freedom represented by the field 
Qic to end up with an effective action in terms of the 
electromagnetic potentials and the order parameter only. 
To this end, one needs a parametrization of the Q/c 
field which explicitly resolves the nonlinear constraint 
([T8|. Following Refs. [23ll23 |. we adopt the exponential 
parametrization 



k(r)=e-^W/2Ae^W/2, 



(23) 



where the matrix multiplication in the time space is im- 
plicitly assumed. The matrix A represents the normal 
metal saddle point (hereafter we work at T > T^) in the 
absence of external fields. 



A,,. 







where — g^® and 







® To 



(24) 



(25) 



The function Ft^t' is the Fourier transform of the equi- 
librium distribution function ~ tanh(£/2T), i.e.. 



iT 



sinh(7rTi) 



t>l/T I 

2T 



S'{t). 



(26) 



S-iQ, A,A,$] ^Sa + S^ + S„ 



(20a) 



The last expression is an approximation applicable for 
slowly varying external fields. Note that choosing the 
parametrization in the form [Eq. (j23p ] does not imply 
that the electrons are in the state of the global thermal 
equilibrium. Indeed, the actual distribution function is 
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given by e'^'='''('-'*)Ft_t/e-*^'='''(''^*'), cf. Eq. ((T9]), and in- 
cludes local variations (e.g., chemical potential) due to 
the presence of electromagnetic potentials. The field K, 
is to be chosen (see below) to achieve this goal in an 
optimal way. 

The matrix field Wt.v (r) in Eq. represents fluctua- 
tions of the electronic degrees of freedom and is to be inte- 
grated out. To avoid redundancy of the parametrization, 
one needs to ensure that the matrix W does not commute 
with A. This is achieved by requiring that VI^A-f AW^ = 0. 
This condition is resolved by introducing four real fields 
wj^, (r) , w^, (r) with a = 0, 2 representing diffuson de- 
grees of freedom and two complex fields (r), iDtc (r) 
for Cooperon degrees freedom. (The bar symbol denotes 
an independent field, not a complex conjugation.) These 
fields are built into the matrix— 



— W*T^ WqTq + WzTz 
WqTq + WzTz WT+ - W*T- 



(27) 

where the asterisk stays for complex conjugation and the 
matrix U is defined in Eq. (|25p . 



B. Diffuson modes, A3— gauge, and normal action 

In this subsection, we shall disregard the fluctuations of 
the order parameter A. Since we are not interested in the 
weak-localization effects, we can disregard the Cooperon 
degrees of freedom wu' (r) , Wu' (r) in the matrix (P7|) as 
well. We then substitute the matrix W, Eq. (|?7)) . written 
in terms of the diffuson fields wf^/ (r) , wf^, (r) (a = 0, z) 
into the sigma-model action (|20cp and expand it to the 
linear order in the diffuson fields. We focus first on the 
z-components ^w^ . Demanding that the terms linear 
in w^,iju^ vanish, one obtains the condition^^ (for details 
see Appendix [X)) 



A] + D{EVA^ - ASVAkA) = 0. 



(28) 



This matrix equation may be resolved by a proper choice 
of the gauge doublet K- ~ (A^"^', /C^), thus fixing the /C- 
gauge. Employing Eqs. (|22ap and ((24)) . one may rewrite 
Eq. ([^5]) as an explicit gauge-fixing condition 



$/c(r,w) = 



1 



-2B. 
-1 



DdivAK:ir,u;) 



with 



= coth(w/2r) 2T/LU 



(29) 



(30) 



being the equilibrium bosonic distribution function. In 
the absence of the quantum components of the fields 
(used to generate observables), Eq. (^5]) is reduced to the 
gauge condition ([7]) written for the classical field compo- 
nents. It is important, however, to fix the gauge for both 
quantum and classical components. 

Equation completes the task of finding the gauge 
field and combined with Eqs. (fTO|) and (|24p provides the 



approximate saddle point, which is determined for any 
given realization of the fields $(r, t) and A{r,t). This 
general scheme guarantees that in the expansion over W- 
fluctuations, terms such as Tr[<I>;cW^] and Tr[Aj<;W^] do 
not appear in the action. 

This procedure does not completely eliminate terms 
linear in the diffuson generators. Indeed, contributions 
of the form Tr[AA^cAW^Ax;] come from the diamagnetic 
term of the cr- model action (j20cp i^ Such terms are lin- 
ear in the diffuson flelds , and quadratic in the elec- 
tromagnetic potentials. Integrating out diffusons w'^,w'^ 
in the Gaussian approximation yields a nonlocal vertex 
quartic in the electromagnetic potentials. It is exactly 
this quartic vertex which is responsible for Altshuler- 
Aronov correction to the conductivity of normal metals^^ 
(for details see Ref. [23| ) • Since Altshuler- Aronov correc- 
tions do not exhibit a singular temperature dependence 
in the vicinity of Tc, we shall ignore these terms hereafter. 
It is an interesting and open question to investigate other 
possible implications of these nonlinear terms. 

Once the terms linear in diffuson generators are elimi- 
nated by a choice of the proper gauge, one may substitute 
the metallic saddle point = A into the sigma-model 
action [Eq. (PO)) ] to obtain the effective action in terms of 
the electromagnetic potentials (for details see Appendix 
\K\ . Such a procedure neglects nonlinear interactions of 
the diffuson modes and thus amounts to disregarding the 
Anderson localization effects. The resulting action takes 
the form 



S 



N 



A^Tz^Ak 



where the operator 7]j is defined as 

Jd = 



0_^ -dtD 
Ddt AiTD 



(31) 



(32) 



The arrows ^ on top of the time derivative imply that 
the differentiation is performed to the left/right, respec- 
tively. Employing gauge fixing condition [Eq. (|29p]. the 
action may be rewritten as 



S 



N 







AiT 



(33) 

This action may now be employed to determine (fiuc- 
tuating) charge and current densities. To this end, 
one needs to introduce an auxiliary vector Hubbard- 
Stratonovich field (r, t) , to decouple the term quadratic 
in quantum component of the vector potential Ajc 



e^uDT Tr[A^] ^ _ 



-(4e^iy_DT)"^Tr^?-2iTrAJ.^j 

(34) 



The resulting action is now linear in both and A^ 
fields, allowing to define the charge and current densities 
as 



, , 1(55[$,A] 



2 6^i{r,t) 



(35a) 
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J(l-,i) = 77 ■ 



2 (5A9(r,t) 



(35b) 



It is important to note that the difFcrcntiation here has 
to be performed over the bare electromagnetic potentials, 
while the action is written in terms of the gauged ones. 
The connection between those {<!>, A} ^ {^!c,A.ic] is 
provided by the functional /C[$, A], which is implicit in 
Eq. (P^)) and in the explicit form is presented in the Ap- 
pendix El A simple algebra then leads to a set of the 
continuity equation (|12p and the expression for the nor- 
mal current density 

j(r,t)=i?[e2z.E(r,t)-Vp(r,t)] -t-Cj(r,t). (36) 

The Hubbard-Stratonovich field ^j(r, t) has a meaning of 
the Gaussian Langevin noise sourco22, with the correla- 
tion function given by (cf. Eq. (|34p ) 



(C?(r,t)^f (r',i')) - 2Te^vD5^^p5t-t'5r-r 



(37) 



Notice that because of assumed local equilibrium of elec- 
tronic degrees of freedom, Eqs. ((55)) and ([57)1 do not lead 



J 



to any excess noise beyond the one prescribed by equilib- 
rium FDT. This is not the case for fluctuating supercon- 
ductors. Indeed, the order parameter may be driven out 
of equilibrium, while the electrons are still in the state of 
the local equilibrium. 



C. Cooperon modes and superconducting action 

Having taking care of the diffuson modes with the help 
of the /C-gauge, we turn now to the fluctuations of the 
Cooperon modes w,w. The latter are induced by the 
fluctuating order parameter A(r, f). To eliminate the 
Cooperon degrees of freedom of the electronic system, we 
substitute parametrization (|27| into the action (|20cp and 
expand it to the second order in w,w. Once again, ne- 
glecting the higher order terms in the expansion amounts 
to disregard the localization effects. As a result, one ob- 
tains the following quadratic action (for details see ap- 
pendix |B]): 



(38) 



here we have introduced vector 2nft'(r) — [wtf (r), (r)]-^, defined in the two dimensional space of the complex 
Cooperon fields and the vector (r) = [Vu'{r), — Vtt'(r)]^, with the elements 



Vw{r) = <5*-,.A^(r,t) + ^6'^_^,AUr,t'), Vwir) = <5,_M^(r,i) - A^(r, i). 



(39) 



The factor iS'(t — t')/(2T) multiplying the quantum component of the order parameter is nothing but the long time 
approximation for the fermionic distribution function Ft-f, Eq. (|26p . This approximation is adopted throughout 
the subsequent calculations. The electromagnetic field-dependent Cooperon matrix propagator has the following 
structure: 







■ Jitv Mw ' 


. cy_ 


+ 


_ Mw y^w _ 



with the matrix elements 

C;,} - -dt + dr + ie[^t(v, t) - $^(r, t')] - [V - zeA^(r, t) - ieA^(r, t')] ' , 

C-} = dt- dr - ze[<i>^(r, t) ~ $^(r, t')] -D[^- zeA^(r, t) - zeA^(r, t')] ' , 



2eD 



■ div A«.(r, t) + A«.(r, t) [V - 2ieA^(r, t)] 



(40) 

(41a) 
(41b) 

(41c) 



Mw 



-2e^D 



A^(r,t) 



2T 



5tA^(r,t) 



KM. 



(41d) 



We show below that the diagonal elements of £ ^ op- TDGL theory in the form derived by GEi^ More pre- 
erator are responsible for the conventional part of the 
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ciscly, the terms and yield the TDGL equation 
for the order parameter, while the additional diagonal 
terms !N" (proportional to the quantum component of the 
vector potential) lead to the superconductive part of the 
current. Interestingly, this is not the entire story yet! 
Indeed, the operator contain also the off-diagonal 
elements M, which induce cross correlations between w 
and w Coopcron generators (i.e., they induce correla- 
tions between rotations of retarded and advanced sectors 
of the Q-matrix in the Keldysh space). It is very dif- 
ficult (if not impossible) to capture such terms within 
the analytical continuation technique of GE. It is exactly 
these off-diagonal terms which are responsible for MT 
contribution to the average current. We shall derive the 
corresponding part of the effective TDGL action, which 
allows to include MT effect in the higher moments of the 
current as well. 

The next step is conceptually simple. It involves the 
Gaussian integration over the Coopcron degrees of free- 
dom (i.e., vector 211) in Eq. ((38)) . The result may be 

schematically represented as S" ~ Tr[y^^£[Ayc] V], which 
is an action quadratic in the order parameter A. One 
needs thus a way to invert the operator given by 
Eqs. (|40)) and ([4T|) . To this end, we notice that the sec- 
ond term in the rhs of Eq. (|40p contains the quantum 
component of the vector potential and thus may be re- 
garded as small. Taking advantage of this fact, we first 
find the saddle point of the action ((38)) without the last 
term in the rhs of Eq. (|40)) and then substitute this ap- 
proximate saddle point into the J4 and M terms. Taking 
variation with respect to w and w, one finds for the ap- 
proximate saddle point 



w. 



.(r)^- dr'dr/C-l^,V_^.^{r'), (42a) 



Wr+i,r-^{r) = 2* |drWC^;;;,,,F^v,,_:^(r') (42b) 

and similar expressions for conjugated Cooperon gen- 
erators w* and w* . The retarded (i] > rj') and ad- 
vanced (?7 < 7]') Cooperon propagators are determined 
by Eqs. (fT5)). 

We now substitute the saddle point Eq. (|42)) back 
into the action (|55)) . to obtain the effective action of a 
fluctuating superconductor (details are presented in Ap- 
pendix |B|) 

5cff[AK;,$K:,AK:] = Sgl[Ak:, A^:] + Ssc[^k, Ak:] 

+ Smt [Ak ,A^] + Sn[^k,A^]- (43) 

The last term here Sn is the normal action ()3H) , or cquiv- 
alently (|55)) , originating from the diffuson degrees of free- 
dom. The other three terms originate from the Cooperon 
action (|38)) in the way outlined above. Specifically, the 
Ginzburg-Landau action Sgl comes from the C~^, 
terms in the action, the supercurrent action Ssc from 
the 3\f terms and Maki-Thompson action Smt originates 




A(r.t) 



FIG. 1: (Color online) Diagrammatic representation of the 
effective action, Eq. (|43p . (a) Conventional Ginzbirg-Landau 
functional Sgl- (b) Anomalous GE coupling between the 
scalar potential and the order parameter, (c) Gradient part 
of supercurrent action Ssc- (d) Diamagnetic component of 
the supercurrent. (e) Nonlocal MT term (here t — {t + t') /2 
and there are two possibilities: one vector potential classical 
and another quantum, which is part of the current, or two 
quantum, which is FDT counterpart). 



from the off-diagonal M terms. The diagrammatic rep- 
resentation of these terms is given in Fig. [TJ 

The TDGL part of the action SGL[^K,A'fi] has a local 
form standard for the Keldysh formalism. It comes from 
Tr[\7tcy] and Tr[ytcy] terms. Fig. (Ha) 



GL 



■nv 



(44) 



where the fluctuations propagator L has a typical bosonic 
form in the Keldysh space 



L-/ 



(45) 



Here, retarded and/or advanced components of the fluc- 
tuation propagator are given by 



Tdt 



Tal + D[\7~2^eA^^ir,t)f 



(46) 



while Keldysh component of the propagator satisfy the 
FDT in cquihbrium Lj} = (L^^ - L^^) AiT if 
ijj <^ T. Note that the scalar potential although 
present in the action ()38p through the operators ()4Iap 
and (|4Ib|) . does not show up in the Ginzburg-Landau 
action (j^ . This happens because upon substitution of 
the Cooperon generators by their saddle point values (|42)) 
the terms $^(r,t) and $^(r,t') in Eqs. (jila)) and (|4Tb)) 
cancel each other. [To be precise there is a small residual 
term ^ dt^'^/T, which we do not keep, since it exceeds 
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the accuracy of our calculations. For the same reason, 
terms with are not kept in the DSf operator Eq. (|41cp ]. 
As a result, the effective TDGL action depends only on 
the vector potential, but not on the scalar potential, if 
written in the /C-gauge. In any other gauge, there is a 
linear coupling between the scalar potential and the z- 
diffuson mode ~ i^^'^'w^ + cf. Eq. (gS]). Taken 

together with the terms AiDw^ and Aww^ (see next sub- 
section) and being averaged over the diffuson fluctua- 
tions, these terms lead to a nonlocal coupling between 
the scalar potential 4> and the order parameter |Ap, see 
Fig. [IJb). This is the anomalous GE term. Thanks to 
the condition (|28p the anomalous term is absent in the K,- 
gauge, making this gauge especially convenient to work 
in. 

The supcrcurrent action Ssc [^K ■, ^k] comes from the 
diagonal terms in Eq. (j38p Tr[w?sfi(j] and Tr['(i;]sr'!Zi], where 
Cooperon generators w and w are given by Eq. (j42p . It 
is also local and given by 



Ssc 



TievD 
2T 



Tr{Apm [A;^^'(V - 2*eA^)A^] } 



(47) 

The gradient part of the supcrcurrent originates from 
\ div A^ -I- A^V terms in the rhs of Eq. \Alc\ . while the 
diamagnetic current originates from 2ieA^A^ term, see 
Figa.IHc) andlHd). 

Time locality of the Ginzburg-Landau and the super- 
current actions stems from their diagonal nature. Indeed, 
they both involve products such as wt^t'Wf j- According 
to Eqs. dSll) and (O, wt^t' ~ 0{t - t'), where e{t) is 
the step function smeared at the scale 1/T. Therefore, 
the time variables in Wt^t'Wf^t are compatible only in the 
narrow vicinity oi t = t', i.e., \t — t'\ < hence the 

time locality. This argument does not apply to the off- 
diagonal MT term. Indeed, the latter involves product of 
retarded and advanced generators Wt^t'Wt',t whose time 
variables are compatible for any t > t' . 

The Maki-Thompson action S'7\/t[A^, A^;], coming 
from the off-diagonal blocks of operator, has essen- 
tially time nonlocal form as explained above, Fig.[T]^c), 



Smt 



e^v Tr 



A;^(r,<)T5i,(t,t')AK(r,t') 



(48) 



where the operator 7sD{t,t') is given by [cf. Eq. (|32p ] 



"^SD — 







MT 
r.t.t' 



(49) 



where the (5Z)^^"'"[A^] functional is given by Eq. 
Note that the MT action has exactly the same structure 
as the second term in the normal action ([5T|) . It there- 
fore amounts to the time nonlocal renormalization of the 
normal diffuson constant D5t-t' DSt-f + 5D^^^^,. 

Finally, we comment on the so called density of states 
(DOS) contributions!^ They originate from the sublead- 
ing terms (in characteristic frequency over temperature) 
in the diagonal operator J-i [not written explicitly in 



Eq. (|41cp ]. see Appendix IB 31 for details. Accounting 
for them leads to a local renormalization of the den- 
sity of states prefactor in the normal action, Eq. (|3ip . 
or Eq. 



1 - 



7C(3) I 



(r,t)| 



(50) 



where C,{x) is the Riemann zeta function. This is a small 
effect in the regime we are working in. 



D. Nonlinear terms in the Ginzburg-Landau 
functional 

The nonlinear in A terms of the TDGL equation are 
not very significant at T > T^- Nevertheless, we shall 
discuss them here for completeness. There are several 
ways "|A|^" terms appear in the effective TDGL action. 
We shall keep track of A*''A'^'|A^'|^ terms which directly 
contribute to TDGL equation for A"^', discussing other 
combinations only briefly. The most important way such 
terms appear is through the third order expansion of the 
Tr[Aj(cQK:] term of the cr-model action (|20cp in powers of 
W . Keeping only the Cooperon generators and employ- 
ing Tr[T-|-r^T:^T±T:^] = =F 1 5 0^0 obtains 



Snl = Tr [A^ {w*ww* — w*ww*) + c.c. 



(51) 



Similar terms coming with A^ component of the order 
parameter eventually cancel out between w and w con- 
tributions and thus are omitted in Eq. ([5T|) . Next, one 
substitutes the saddle point value of the Cooperon gen- 
erators Eq. ([42]) into the action (jSTjl . and perform traces 
over the time indices. In doing so, one should keep only 
the first power of the quantum component of the order 
parameter [coming from Eq. p9p ]. The diagrammatic 
representation of the corresponding terms is shown in 
Fig. mja). After straightforward algebra (see Appendix 
[C]for details), one finds 



S. 



NL 



7C(3)^ 
87r2T2 



Tr[A^''(r,t)A^(r,t)|A^(r,t)p+c.c.: 



(52) 

This term is to be added to the retarded and advanced 
(but not Keldysh) parts of the Ginzburg-Landau action 
P^ . It leads to the standard nonlinear term of the 
Ginzburg-Landau equationf22i2^ which is (a) local, (b) 
disorder independent, in agreement with Anderson theo- 
rem. 

Interestingly, this is not the only way the nonlinear 
terms appear in the effective action. Let us mention two 
other venues, (i) One may expand the D Tr[V(5;c]^ term 



of the nonlinear a model action (|20cp to the fourth order 
in the Cooperon fields, generating the so-called Hikami 
box, see Fig. [UJb) . One then substitutes the saddle point 
value of the Cooperon generators (|^^ in such a term 
to obtain a contribution quartic in the order parameter. 
This term leads to a local renormalization of the diffusion 
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K'"{r,ty 




h''(r,t) 



A'i{r,t) ^'''{r,t) 




(d) 



FIG. 2: (Color online) Nonlinear contributions to Ginzburg- 
Landau functional, (a) Diagram for the local nonlinear term 
of the action Ea. (|52|) . (b) Renormalization of the diffusion 
constant in the TDGL equation - shaded square is the Hikami 
box. [(c) and (d)] terms involving diffuson channel, leading 
to a nonlocal renormalization of the dynamic part in TDGL 
equation. 



constant in TDGL equation SDr,t ~ £)|A='(r, <)| V^^- 
There is no MT nonlocal renormalization of the diffusion 
constant in the superconductive part of the action (as op- 
posed to the normal one where both MT and DOS renor- 
malizations take place), (ii) There is yet another source 
of nonlinear terms (we are not aware if it had been dis- 
cussed previously in the literature) . It originates as a re- 
sult of mixing between Cooperon and diffuson channels. 
To see it, one expands Tr[A;c'3K;] term of the c-model 
action (|20cp to the second order in W. This way one 
generates interaction vertices of the following structure 
Tt[Aicw* (w" ± w")] + CO. with a = 0, z and correspond- 
ing terms with w. We then perform Gaussian integration 
over the diffuson fields , w" to obtain a nonlocal vertex 
TT[w;^,AK:ir,t'){wf,^wf^,)A*^{r',i')wt,f], where is 
the diffuson propagator, Eq. (jlip . There is a similar ver- 
tex with w*,w generators. Such a nonlocal vertex is ef- 
fectively a renormalization of the diagonal part of the 
operator in Eq. ([38|) . It is important to stress that 
the diffuson admixture does not generate the off-diagonal 
terms in operator, thus not affecting directly the MT 
channel. We then substitute the saddle point values of 
the Cooperon generators, Eq. (H^ . in this nonlocal ver- 
tex and find 



S 



= S^Tr [A*^^r,t)A'fl{r,t)Vl:^:Al^\r',t')d,A^^{r',t')+c.c.] 



(53) 



r 



This term is formally of the same order of magnitude 
as the conventional one, Eq. ([5^ . Indeed, the diffuson 

propagator Vl'^, is not cut by the temperature [unlike 
the Cooperon propagators in Fig. [lb] . 



E. Equations of motion and Coulomb interactions 

To derive the stochastic equations of motion presented 
in Sec. [ill one needs to get rid of terms quadratic 
in quantum components of the fields: in Sql and 
Smt- This is achieved with the Hubbard- 



in Sn 



Stratonovich transformation similar to Eq. 
and 



for 



(54) 

for A^. As a result, the effective action (^5)) acquires the 
form linear in quantum components of the fields. Integra- 
tion over the latter leads to the functional delta functions 
imposing the stochastic equations of motion. This way 
the TDGL equation ([5]) , which we present here including 
the nonlinear terms 



J 



dt+Tal D{V 2*eA^(r,t))2 + ^|A^(r,t)p -^j dt'dv' Al{v' X)drA^{v' ,t')Vl% 



AK(r,t) =CA(r,t) (55) 



with h = 7C(3)/7r^, the continuity equation (|12p . and noise sources may be directly read out from the Hubbard- 

expression for the current (jlSp arc obtained for the clas- Stratonovich procedure and are given by Eqs. ([5]) and 

sical components of the fields (we have omitted the sub- (fTB)) . 
script cl for brevity). The correlators of the Gaussian 
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We discuss briefly the role of the last nonlinear term 
in the Ihs of TDGL equation ([55]) . On the mean-field 
level, i.e.; being averaged over the fluctuations of the or- 
der parameter (A^(r, i)A;^='(r', t')) (x Lxir - r' ,t - t'), 
this term leads to a renormalization of the coefficient in 
front of the time derivative IHl), 



dtAK 



4-d 



dtA 



(56) 



where = VD/Tc and ^{T) = ^D/{T-T,). The 
dimensionality dependent coefficient Cd appears as the 
result of the convolution between the diffuson and the 
Keldysh components of the fluctuations propagator and 
reads as 



Cd 



7C(3) 
87r3 



1 



[2ttY (fc2 + l)(fc2 + i/2)' 



(57) 



here is ci = 0.044, C2 = 0.012, and cg = 0.005. Note that 
disorder-dependent renormalization of the dtAjc term 
does not violate Anderson theorem. The down renormal- 
ization of the coefficient in front of the first time deriva- 
tive is a precursor of the oscillatory Carlson-Goldman^i 
modes, appearing below Tc- 

Finally, the TDGL equation ([8]), written in the K,- 
gauge may be transformed to an arbitrary gauge by the 
substitution A^; = Ae"^*'^'^. Such a substitution brings 
the vector potential A to the covariant spatial derivative, 
while the the time derivative acquires the anomalous GE 
term dt ^ dt — 2iedtlC. This way TDGL equation ^ is 
obtained. 

The set of equations is simplified in the limit of the 
strong Coulomb interactionsi^ The latter impose the con- 
dition of local instantaneous charge neutrality p{r, i) = 
and therefore div j = 0. Applying this condition to the 
expression for the current p3p and neglecting for sim- 
plicity the MT and supcrcurrent contributions, one finds 
e^vDdvv E = div ^j. As a result, the longitudinal com- 
ponent of the electric filed E|| is a fluctuating Gaussian 
field. Employing Eqs. PH)) and pB]) . one may translate 
it to the Gaussian correlator for the longitudinal com- 
ponent of the vector A«,Q{r,t) ^ A.fJ^{t) + A]|''"''*(r, i). 



where A|j^*(t) is an externally applied divergenceless field 
and the fluctuation component has the correlator 



fluct 



(-q, -a')A(''"=*(q,Lj 



1 



2T 



e^vD (L>q2)2 



Exactly the same expression may be, of course, directly 
read out from the normal action ([33|) (the Maxwell part of 
the action is absent in the the strong Coulomb limit). To 
this end, one needs to perform the Gaussian integration 
over the component leading to the quadratic action 
for the A^ with the correlator given by Eq. ([58|) . It is 
this fluctuating vector potential which is responsible for 
dephasing of the Coopcron propagatorsi^ 



F. Effective action versus diagrammatic technique 

Given that there is an existing microscopic formalism 
for the Aslamazov-Larkin, Maki-Thompson, and density 
of states diagrams, it is important that the results from 
the Keldysh effective action, formulated in the previous 
sections, be compared with well established results for 
the corrections to the conductivity. This comparison is 
the necessary check for the validity of our approach. 

We start from the density of states contribution to the 
conductivity. For that purpose, one uses Eq. (j50p and 
writes conductivity correction in the form 



47r2r2 



(|A^(r,t)|), 



The averaging over the order parameter fluctuations is 
done easily in the momentum space (|A^(r,t)|)^ = 
Lif (q, cj). Using then explicit form of the 
Keldysh component for the fluctuations propagator ((45|). 
one finds 



Sa 



DOS = 



28C(3) 



+ CXD 



•GL 



)-fa;2 



(60) 

Further analysis of the formula ([60|) depends essentially 
on the system effective dimensionality. As an example, 
let us concentrate on the quasi-two-dimensional geometry 
- metal film with the thickness d, which is much smaller 
than the superconductive coherence length d <^ ^{T). In 
this case, momentum sum can be written as the integral 
according to the substitution Ylq ~* ^ / Perform- 
ing remaining integrations, with the logarithmic accuracy 
one finds 



Sa 



DOS 



7C(3)e2 



In 



T-Tr 



(61) 



Deriving Eq. (|6ip . the momentum integration was cut at 
the upper limit Dq^^^ ~ Tc- Recall that effective action 
was derived under the constraint Dq2 ^ cj <c T^, thus 
such a regularization is self-consistent. 

We proceed with the Maki-Thompson correction to the 
conductivity. In this case, one starts from the formula 
Sctmt = {^D^'^) uses expficit form of the SD^'J^ 
given by Eq. (fH|) . and rewrites average over A in the 
momentum space. This way one obtains 



5(Jmt = — e^D J duj 



l~l2 I ,,,21 



(62) 

Again in the case of the two dimensional geometry, after 
momentum and frequency integrations, Eq, (|62p reduces 
to 



S<7mt — —, 



8d \ T-TaJ"\TGL 



In 



(63) 
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where infrared divergency in the momentum integra- 
tion was cut here by hand introducing dephasing time 
ZJq^ijj r^"'^. This spurious divergency is very well 
known feature of the Maki-Thompson diagram. It was 
regularized by Thompson introducing magnetic impuri- 
ties, and in that case dephasing time is nothing else 
but the spin flip time Ts- 

Finally, we summarize with few additional remarks. 
Equations ([60|) and (|62p can be recovered from the tra- 
ditional Matsubara diagrammatic techniques after one 
expands all fluctuation propagators at small frequencies 
and momenta, integrates fast fermionic energies, and 
keeps only contribution from zero Matsubara frequency 
[in our language, latter condition strictly speaking cor- 
responds to the long time approximation for the distri- 
bution function Eq. ipS]) ]. To this extent, effective action 
approach contains the most divergent temperature part 
of the conductivity corrections, thus allows to reproduce 
known results. 



IV. DISCUSSION 



supported by the A. P. Sloan foundation. 



APPENDIX A: DIFFUSON EXPANSION AND 
NORMAL ACTION 



We first focus on the part of the action ((20l) which is 
linear in the electromagnetic potentials. There are two 
such terms: originating from the trace TtI^/cQk] and 
5a coming from the covariant derivative (PT|) . We then 
expand Qk: matrix to the linear order in deviations from 



the saddle point Qk 



and use the 



cyclic property of the trace operation to obtain 5*$ = 
_-|iiTr([$K,A]Ty) and 5a = {AVW[EAic, A]) 

which after the integration by parts translates into 
5a = -^^Tt{EVAic~'AEVAk:A)W. Requiring 
that terms linear in variation W (and linear in poten- 
tials) vanish, one arrives at Eq 

Expanding action pO]) to the second order in W and 
keeping track of the diffusons only, one finds quadratic 
action of the diffuson degrees of freedom 



We have presented a systematic way to integrate out 
the fermionic degrees of freedom in a fiuctuating super- 
conductor. The underlying assumptions for this proce- 
dure are: (i) spatial and temporal scales of all bosonic 
fields are slow in comparison with and 1/Tc corre- 
spondingly (but not necessarily slow in comparison with 
C(T) = ^oy^Tc/(T - Tc) and tgl)- (ii) The external elec- 
tromagnetic fields are sufficiently small (see Sec. Ufor the 
details), such that the fermionic degrees of freedom are 
in local equilibrium. The result is the dynamic Keldysh 
action written in terms of the fiuctuating order parame- 
ter and electromagnetic potentials. This action naturally 
incorporates (time nonlocal) MT terms as well as anoma- 
lous GE terms, effectively closing the discussion whether 
or not TDGL theory includes the MT effect. We have 
also uncovered certain nonlocal nonlinear terms of TDGL 
equation (passed previously unnoticed, to the best of our 
knowledge). The nonlinear coupling of the electromag- 
netic fields, leading to the Altshuler-Aronov effect, may 
be also directly incorporated into the scheme. 

Although we did not evaluate any physical observable 
here, the derived action opens the way to describe a num- 
ber of phenomena. To name a few, we mention, e.g., 
nonequilibrium current noise in proximity to the criti- 
cal temperature, especially the MT contribution to noise, 
which would be very difficult to calculate by any other 
mean. Another possible application is evaluation of the 
MT dephasing time. Extension of the theory for T < Tc 
to describe, e.g., the collective modes^i and analysis of 
the dynamical regimes far from the equilibrium, are yet 
another fascinating directions. 

We are grateful to L. Glazman, B. Ivlev, M. Khodas, 
I. Lerner and A. Varlamov for useful discussions. This 
work is supported by the NSF Grant Nos. DMR 02- 
37296, DMR 04-39026, and DMR 0405212. A.K. is also 



iSd[w'>, w^] = — Tr {w'^,, + dt + , 

(Al) 

where a = 0, z. This action leads to the following prop- 
agator: 



..a t„\„-a ( „^\ _ ^ (2"') ^ei-e4^e2-e3 (A2) 



«,e.(qX.e.(-q)) = - 



Substituting the saddle point A in the cr-model action 
pi)) , one finds 



Sn = — Tr{$KT$K} 



(A3) 



-Tr {SAk;ASAk;A-SAkSAk;} 



The traces are evaluated using the following matrix iden- 
tity: 

+ OC 

y d6Tr[T„T;3-T„A,^T;3A,_] =8c^(flji)„;3, (A4) 



where 



n: 



i_^o -1 

1 



(A5) 



and Ti = (To <8) Tz, T2 = (g) r^, e± = e ± bjjl. The 
identity is based on the relation between the bosonic and 
fermionic distribution functions: 



-t-00 

/ 



F,J = 2c^, 



(A6a) 
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+CXD 

J de(l-F,^F,J = 2wB, 



(A6b) 



Combining Eqs. (jA3[) - (|A6[) . we arrive at normal metal 
action (|3ip . Finally, to find the macroscopic equations 
and ([55]) . one needs to express gauged electromag- 
netic potentials ^jc, in terms of bare ones $, A. Using 
Eq. ([29]) . one may find the relation between those 

A^{q,u;) = A^(q,c^) - iqV^flZ'[$ + c^(qA||)/q2], 

(A7a) 

$,ciq,uj) = Dq'vy^A^ + Lu{qA«)/q% (A7b) 



where 



2^2 



-iuj + Dq^ 







(A8) 



APPENDIX B: COOPERON EXPANSION AND 
EFFECTIVE ACTION 

Carrying out the Cooperon expansion, it is convenient 
to distinguish several contributions into the cr-model ac- 
tion (1201): 

5,[M^,A^,$^,A^] = S'',[W]+Sl[W,K^] (Bl) 
+ Sl[W,A^]+Si[W,A^]. 

The 5'^[VF] part corresponds to the free Cooperons 
which are uncoupled from both /S.]c{Y,t) and Aic{Y,t). 
This contribution arises after one expands trace of 
the gradient term Tr(VQ)2 = Tr(AVW^AVM^) = 
TrCWV^W) and trace of the time derivative term 
Til^dtQ] = ^TT[EdtAW^] = ^TT[EdtUAoWWU]. Mul- 
tiplying matrices, tracing them over Keldysh-Nambu 
space, one finds 

= — Tr{<,,[-i?V2 -dt + dt']wt,t + iw^ w)}. 

(B2) 

The ^^[M^, A/c] part corresponds to the coupling term 
between Cooperons and the order parameter. This 
contribution arises from the trace Tr|Ax;(5A:] after one 
expands Qic to the first order in W: Tt[A!cQk] — 
Tr[AK;AW^] Tt[UAk:UAoW]. After evaluation of 
traces, which is done with the help of the identity 
Tin [t±TzT^] = T1; one finds 



St[W,A,c]^^TT 



(B3) 



where we have used vectors 211 and V in the notations of 
Eq. (IMl). 

The c and d parts of the action are the terms which 
provide the interaction vertices between the cooperons 
W and the vector potential A}q. The ^^[VF, A^] part is 
linear in the vector potential and arises from the square 



of the covariant derivative, keeping terms linear in A^; 
contribution 

SI = —-^y{VQkX^AkMk\} = 



^^Tr{AVM^[SA;c,AM^]} = --^Tt{UEAk.U[VW,W]}. 
After a straightforward algebra, one finds 

S-,[W, A,c] —TT{Au"[wu'^y^t't" - Vw^t^wrt.] 

(B4) 

+At"K*t'VWt't" - VWtt,Wt't"] - c.c.} , 

where we have introduced notations 

AtAv) = St^t'A^^{r,t) + ^6[,_tAl{r,t) , (B5a) 



Au'{r)^6t^t'A^^{r,t) - —S[,_,Al,{v,t') . (B5b) 

The quadratic in vector potential part of the action 
comes from the diamagnctic term, which has the form 

S^IQk] = -^^^Tv{EA^QrcEA,cQK}- One possibil- 
ity is to expand one of the Qk: matrices up to the second 
order in W, while leaving the other to be A, the other is 
to expand both of them to the first order 



S'AW,Ajc] 
where 



me vD 



Tr [ttWW + twtw] , (B6) 



(B7) 



F,Al A^(1-F,F,0 + A^(F,-F,, 



\cl 



■ F^i A^ 



and all vector potential have e — e' as an argument. Com- 
bining it with all the contributions, given by Eas. (|B2p - 
(|B4[) , the full action may be conveniently presented using 
matrix notations as Eq. ([55)1 . 

We turn now to the derivation of the effective ac- 
tion, Eq. ([431). To this end, we transform Eq. (|42]) into 
the energy-momentum representation, then using explicit 
form of the V vector, given by the formula ((5^ . we find 
for the Cooperon generators 



Wee'(q) = -2i 



,.A^(q,£-£') + J^eAHq,£-e') 



D{q-2eAf) -i(£ + e') 



D {q-2eAf) +i{e + e') 



a) 



(B8b) 
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Note that the scale for the energy center of mass e = 
{e + is set by the temperature e ^ T. Then m 

most of the cases Dc^ ~ ~ T — Tc can be ignored as 
compared to e (the exception is MT term), thus one may 
write instead of (|B8p the approximations 



'(q) 



e + e' 



zO 



(B9a) 



wtv (r) - iOit' - i)A^(r, r) - A^(r, T)Y{t - 1'), (BlOb) 
where t = {t + t')/2 and 



f de 


" e " 


/ — tanh 




/ 27r 


12t\ 



e + iO TT 



-arctanh 



-TTT\t\ 



(BU) 



Wee' (q) ~ 2 



A^(q,£-eO-fe'A^(q,e-gO 
e + e' -iO 



(B9b) 



and similar equations for the conjugated fields. After the 
inverse Fourier transform, one finds 

wwir) = -z0(t-i')A^(r,T) + A^(r,T)r(t-t'), (BlOa) 



!• Sgl[Ak;, A^] part of the effective action 

Let us concentrate first on the part of the action which 
corresponds to the diagonal blocks C^^ and C^^. These 
two give identical contributions to the action (j38p . thus 
accounting for an additional factor of 2. We find 



-2-Kiu Tr 



AC/ 



F A'^ 



Z) (q - 2eA'=' 



2ie 



(B12) 



where we have introduced energy integration variables as 2e = e + e', cj = e — e', e± = e ± uj/2, and A;tc|± 



A/c(±q, iw). We point out here that contribution to the iSa[Aic] with two classical fields ~ A^A^'^' is identi- 



cally zero. This is manifestation of the normalization condition within the Keldysh formalism. Adding the term 
— 27ri/Tr | A^|^Ajq_/[_D (q — 2eA^)^ — 2ie]| = (due to energy integration of retarded propagator), we obtain for 

5,[A;c,A^] + 5a[Ak] = 5gl[Ak,A^] [cf. Eq. 



5gl [Ak, A^] =;.Tr 



where we have introduced superconductive fluctuations propagator in the form of the integral 



L 



R{A) 



1 

A 



i?(q-2eA^) 



2ie 



(B13) 



(B14) 



In what follows, we show that the latter can be reduced to the standard form given by equation (j46p . Indeed, changing 
e = e — uj/2, adding and subtracting term at zero frequency and momentum we write 



A F- ■ 

de I 

2e 



de 



Fe 



^(q-2eA^) 



iuj — 2ie 



Fe_ 

2e 



(B15) 



where the logarithmically divergent integral in the above formula was cut in the standard way by the Debye frequency 
LUiy. Introducing dimcnsionless variable x = e/2T, integrating second term in the right hand side of Eq. (jB15[) by 

oo 

parts with the help of identity J da; ln(a;)sech^(a;) = — In where 7 = with C = 0.577 being the Euler constant. 



and using the definition of the metal-superconductor transition temperature Tc 



T 2 



dx 



tanh(a;) 



tanh(a::) 



D(q-2eA^')2-ii. 
4T 



Using series expansion 



tanh(a;) = ■ 



2x 



a;„ = 7r(n+l/2). 



A^i' , we have for Eq. (jBlsp , 
(B16) 

(B17) 
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interchanging order of summation and integration, integrating over x with the help of 



dx 



xdx 



ITT 



D(q-2eAg)2-i 



4T 



D(q-2eAg)^-»i^ 
4T 



(B18) 



r 



and recaUing the definition of the digamma function 



,Pix) = -C-J2 



ri=0 



one finds 



T 



n + x n + 1 



(B19) 



(B20) 

Finally, expanding digamma function, using ij/ {1/2) ~ 
7r^/2, and transforming back to the real space and time 
representation i:)(q-2eA^)^-za; -D(V - 2ieA^)^ + 
dt, one derives Keldysh version of the Ginzburg-Landau 
action in the form (|44p with the fluctuations propagator 
given by Eq. 



Ssc[A'i,AK] and Sa/t[A^,Ak] parts of the 
effective action 



Now we concentrate on the contributions to the effec- 
tive action (|43p . coming from the 3sf terms of the matrix 
(|40p , proportional to the quantum component of the vec- 
tor potential. These contributions translate to the action 
of the form 



-Tr \wt,'}^tt'Wt't + Wff,Jitt'Wt't\ ■ 



(B21) 

One uses explicit form of 'N given by Eq. (|4T|) and makes 
use of approximations (|B10[) . As soon as J^tt' cx: St-f, 
one may integrate over t' using regularization J dt'S{t — 
t')e{ t' - t) = 1/2. Note that deriving Eq. gll from 
Eq. (jB2ip . we kept only classical components of the field 
A in the Cooperons (jBlOp . The quantum components 



generate interaction vertices like tr[A^ VA^'^'] , having 
more then one quantum field, are smaller than Eq. (|47p 
by the parameter 1/Ttql «C 1. Indeed, one sees from 
Eq. (jBSp that A^ comes in the combination with the 
fermionic distribution function which according to the 
approximation brings additional smallness by one 
extra power of temperature in the denominator, which is 
in contrast with the term having A^. 

In the similar fashion, one derives Smt part of the 
action. We start from 

5mt[A^, Aye] = -^Tr [wl^,Mtt'Wt.t + w*u,Mtt'Wt't] ■ 

(B22) 

At this point, we again make use of approximation (jBlOp . 
Observe that in contrast to Eq. (jB2ip . where we had 
product of either two retarded or two advanced Cooperon 
fields, which restricted integration over one of the time 
variables, in the case of MT contribution (|B22p . we end 
up with the product between one retarded and one ad- 
vanced Cooperon and the time integration running over 
the entire range t > t'. Precisely, this difference be- 
tween Eq. (jB2ip and (jB22p makes contribution Ssc to 
be local, while Smt nonlocal. Finally, in each of the 
Cooperon fields w, Eq. one keeps only contribu- 

tion with the classical component of the order parameter. 
The quantum component is again smaller by the factor 
of 1/Ttgl « 1. 



3. Density of states contributions to the 
Ginzburg-Landau action SdosIA^, Ajc] 

There are two ways subleading DOS contributions ap- 
pear in the effective Ginzburg-Landau action. The first 
one, not written explicitly in Eq. (|38p . is 



S. 



DOS 



[W^,Ak,$k,Ak] = ^Tr 



S^po,? 



(B23a) 



A^(r,t)[A^(r,i) 



Ft'. 



(B23b) 



Note that in order to reproduce correctly DOS contributions one cannot use the approximate form of the fermionic 
distribution function. In what follows, we deal with the part of the action (|B23p having one classical and one quantum 
components of the vector potential. The other one, having two quantum fields can be restored using FDT. To this 
end, we substitute Cooperon generators in the form (jBlOp into the action (jB23p . We keep only classical components 
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of A (the quantum one produce insignificant contributions) and account for an additional factor of 2 due to identical 
contributions from w and w Cooperons. Changing time integration variables t — t" = t and t + t" = 2r], one finds 



S. 



DOS 



iire vD Tr 



A^(r, in + r/2)[A^(r, ,7 + r/2) - A^(r, ,7 - T/2)]Fr 



(B24) 



x0(,7 + r/2 - t')e{t' - 77 + t/2)A^'=' (^r, !l±lZ^^ [1 



77 - t/2 - t' 



Note that due to the step functions, integration over t' is restricted to be in the range 77 + r/2 > t' > 77 — r/2. Since 
Fr is a rapidly decreasing function of its argument, the main contribution to the r integral comes from the range 
r ^ 1/T ^ 77. Keeping this in mind, one makes use of the following approximations: A^(r,77 + r/2)[A^(r, 77 + 



r/2) - A^(r,77 - r/2)] « rA^(r, 77)a,A^(r, 7;) and A^^' r 



*cl v+T/2-t' 



^K. y-^ — 2 — j^Ky-^ — 2 — ) ^ 

to integrate over t' explicitly J dt'9{rj + r/2 — t')6{t' — rj + r/2) = t9{t). Using fermionic distribution fimction ([26 
and collecting all factors, we find 



(r,^ 



-r/2-f' 



|A^(r,77)|^, which allows 



00 

5f = ire'vDT Tr [A^ (r, A^ (r, t) | A^ (r, t) p] j - 



sinh(7rrr) 



(B25) 



where we set rj ^ t. Performing remaining integration over r and restoring Sdos ~ -^k-^a; FDT, we arrive at 



(B26) 



r 



with 7d given Eg. ([32|) . The other source of the DOS con- 
tributions is the matrix element 3\f itself, where one has 
to restore fermionic distribution function, relaxing on the 
approximation (|26p . Then in the term oc TT['Wtt''^tt''w^f.,], 
after one uses Eq. (jBSp , we need to keep momentum Dq^ 
dependance of the Cooperon and expand over Dq^/e <C 
1. This produces subleading contribution such as Scff cx 
4^Tr{A«-I?V2A^|A^|2}. As a resuh, the effective ac- 
tion accounting for the density of states suppression may 
be cast exactly into the form of Eq. ([55)1 , where one makes 
the substitution u ^ u + 5v^P^ . 



APPENDIX C: NONLINEAR ACTION Sjvl[Ak] 

In this section, we show how one proceed from Eq. (|5ip 
to Eq. ([5^ . As was pointed out above, one needs to 



keep only contributions having one quantum component 
of the order parameter field. Overall, there are three 
possibilities to do that in each of the Cooperon sectors w 
and w. Moreover, it turns out that contributions coming 
from the w and w are identical, thus accounting for the 
factor of 6. We thus obtain 



Snl ^—TT[^*,^\v,t)wu'{v)w;,^,,{v)wtntiv)+c.c.]. 

(CI) 

We next substitute the approximate form of the 
Cooperon generators, Eq. (jBlOp . into this formula. In 
the case of w* , we keep quantum component of the order 
parameter and in the other w the classical ones. 



Snl 



Tr 



e{t 



mt" -t)Y{t' ~nA*£'{r,t)A'fi (^r, 



t + t' 



t' + 1" 



t + t" 



+ c.c. 



(C2) 



We change now integration variables as i" — < = r and 
t + t" = 2j] and observe that the integration over t is 
restricted to be in the range r] + T/2 > t > 7; — r/2. Recall 
that according to the definition (|Blip . the function Y{t) 
is rapidly falling on the scale r ~ l/T <C 77. Thus, the 



major contribution to the above trace comes from the 
small r. Thus, everywhere except the theta functions, 
one may set t ^ t' ^ t" ~ rj and integrate over t explicitly 



16 



getting t9{t). Finally, using the integral and collecting all factors, one recovers Eq. (|52p 

Jdr T arctanh [e'^^^] = (C3) 
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